Self-assembly and nonlinear dynamics of dimeric colloidal rotors in cholesterics 
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We study by simulation the physics of two colloidal particles in a cholesteric liquid crystal with 
tangential order parameter alignment at the particle surface. The effective force between the pair is 
attractive at short range and favours assembly of colloid dimers at specific orientations relative to the 
local director field. When pulled through the fluid by a constant force along the helical axis, we find 
that such a dimer rotates, either continuously or stepwise with phase-slip events. These cases are 
separated by a sharp dynamical transition and lead respectively to a constant or an ever-increasing 
phase lag between the dimer orientation and the local nematic director. 
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Self-assembly is one of the key aims of present-day nan- 
otechnology. The idea underpinning this concept is that, 
through a careful control of the interactions in a suspen- 
sion of particles, it is possible to drive the spontaneous 
formation of an ordered target structure or pattern, start- 
ing from a disordered initial condition. The process is 
spontaneous as it often entails the minimisation, whether 
global or local, of the free energy of the system. An out- 
standing specific challenge in modern self-assembly is to 
build a structure with a target dynamic feature, such as 
a synthetic microscopic rotor, walker or swimmer 

Colloidal dispersions in liquid crystals offer a useful 
system for the development of self-assembly strategies 
at the microscale. Even in nematics, the simplest liq- 
uid crystalline phase, elastic distortions and topological 
defects mediate a variety of interactions resulting in the 
formation of wires, colloidal crystals, cellular solids, and 
clusters entangled by disclinations [2H4J. This variety of 
self-assembled structures is possible because one can tune 
the liquid crystal mediated interactions by changing, e.g., 
the strength and nature of the liquid crystal ordering at 
the colloidal surface. These alter the local symmetry of 
the director field near the particle and hence qualitatively 
change the effective interparticle forces. 

Here we show that a powerful way to extend the po- 
tential for self-assembly of colloids in liquid crystals is 
to consider cholesterics, or chiral nematics. In cholester- 
ics the direction of the molecular order in the ground 
state spontaneously twists. The spatial modulation of 
the twist is in the micron range, therefore of the order 
of typical colloidal sizes. We have previously shown [H| 
that by varying the ratio between these two fundamental 
length scales, it is possible to change continuously the 
topology of the defects, or disclination lines, surrounding 
a single particle — morphological changes with no direct 
counterpart in nematics. Very recently, Mackay and Den- 
niston [6(] took a step further and studied the interparticle 
elastic force felt in a cholesteric by two colloidal spheres 
with tangential anchoring of the director field at their 
surface. As in the nematic case Q, there is a complex in- 
terplay between repulsive and attractive directions which 



leads to the formation of a dimer or longer chains. Here 
we focus on the simplest case of a dimer\ but progress 
beyond the purely static investigation of [fj to show that 
such a dimer exhibits unexpected and intriguing dynam- 
ical properties. When subjected to an external force (for 
instance gravity) along the cholesteric helix, the dimer ro- 
tates about this axis in a screw-like fashion. Depending 
on the magnitude of the force, the dimer either rotates 
continuously or exhibits phase slippage, alternating peri- 
ods of smooth rotation with static spells in which it trans- 
lates without rotation. This dynamical transition shows 
similar near-critical behaviour to the depinning of driven 
vortices and of charged density waves in superconductors, 
both of which may be studied with the Frenkel-Kontorova 
model for transport in a periodic potential Another 
analogue is provided by the synchronization of coupled 
oscillators described by the Kuramoto model Q . Within 
our liquid crystal context, the phase slippage regime re- 
quires a specific free energy landscape which we discuss. 
This provides potential for the design of self-assembled 
systems with tunable dynamic properties. 

The system we study consists of two spherical colloidal 
particles of radius R moving in a cholesteric liquid crys- 
tal. To describe the thermodynamics of the chiral host, 
we employ a Landau-de Gennes free energy J 7 , whose 
density / may be expressed in terms of a traceless and 
symmetric tensor order parameter Q [lfj and is detailed 
in . Tangential anchoring is modelled by a surface free 
energy, f s = ^W(Q a p — Q^) 2 , where W is the strength 
of anchoring and Q° a a is the preferred order parameter in 
the tangent plane to the local spherical surface [l2j . 

We employ a 3D hybrid lattice Boltzmann (LB) algo- 
rithm [13j to solve the Beris-Edwards equations for Q [10j 
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Here, T is a collective rotational diffusion constant and 
D t is the material derivative for rod- like molecules fioj ]. 
The term in brackets is known as the molecular field, 
which in the absence of flow drives the system towards a 
free energy minimum. The boundary conditions for the 
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order parameter on the colloidal surfaces are given by 
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where z/ 7 is the local normal to the colloid surface. 

The velocity field obeys the continuity and Navier- 
Stokes equation, with a stress tensor generalised to de- 
scribe liquid crystal hydrodynamics [10[ . Within our hy- 
brid scheme, we solve the Navier-Stokes equation via LB, 
and Eq. Q] via finite difference [f| . Colloids are repre- 
sented by the standard method of bounce-back on links, 
which leads to a no-slip boundary condition for the veloc- 
ity field (see 0, E3 for details). Order parameter varia- 
tions create an additional clastic force acting on the par- 
ticle which is computed by integrating the stress tensor 
over the particle surface 0, Q • 

The dynamics is primarily controlled by the Ericksen 
number, Er = jivR/K, where the rotational viscosity 
7i = 2q 2 /T, v is a velocity characteristic of the flow, and 
q is the degree of ordering in the system. In the uniaxial 
case with director n, Q a p = q(n a np — 5 a p/2>). 

In what follows, we quote our results in simulation 
units [III EH- To convert them into physical ones, we 
can specify an elastic constant of 28.6 pN, and a rota- 
tional viscosity of 1 poise. (These values hold for typical 
materials, and a colloidal diameter of 1 ^m.) In this way, 
the simulation units for force, time and velocity can be 
mapped onto 440 pN, 1 fj,s, and 0.07/xm/s respectively. 

We first consider the interparticle elastic potential. 
Two particles are placed centre-centre separation vector 
d apart. Both particles are held fixed for the duration of 
a simulation in which the free energy is minimised. By 
repeating simulations for different d, we map out the ef- 
fective two-body potential as a function of the reduced 
separation (d — 2R) /p in the three coordinate directions, 
and as a function of angle in the x — y plane (Fig. la). 
The potential is markedly anisotropic. While the poten- 
tial in z shows strong repulsion at large separations, there 
is an attraction in the x — y plane. The most favourable 
configuration at small separations is along the director 
field (x-direction in Fig. la). Here we estimate a max- 
imum attractive force of 20 pN. Before this deep mini- 
mum is reached the dimer needs to overcome a repulsion 
(peak at (d — 2R)/p w 0.2), which is largest when the 
disclinations at the opposing particle surfaces join up. In 
this bound state, the colloids share two disclination lines 
which act as a glue between them (right Fig. lb). For 
separation vectors perpendicular to n we find a stable 
minimum and no repulsive barrier (y-direction in Fig. 
la). These results are far from the nematic limit stud- 
ied experimentally in Q and theoretically in [l5j , as well 
as from results obtained in a twisted nematic cell [l6| . 
Most notably, the in-plane potential perpendicular to n 
(here y) was always repulsive in the nematic 15] . The 
preferred configuration (here, along n) is at an angle of 
about 30° in the nematic 0, EH . 
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FIG. 1: [Color online.] (a) Interparticle elastic potential as a 
function of reduced separation (d — 2R)/p, along the helical 
axis (z) and in the x— and y— directions. The inset shows the 
angular dependence in the x — y plane at (d — 2R)/p « 0.05. 
Units are in ksT w 4 pN nm. (b) Snapshots of the dimer 
showing the disclination lines at (d — 2R)/p w 0.05 when d is 
along y and along x. The Cartesian axes, the separation d, 
and the far-field nematic director, n are shown, (c) Definition 
of the angles and 4> used in the text. For error analysis see 
Supplementary Material [TH |. 



These results confirm and extend those of on the 
energetics of dimer formation in cholesterics. Our key 
focus in the present work is dynamics Wc place two par- 
ticles initially near the weaker minimum of the potential 
in the direction perpendicular to n (left Fig. lb). This is 
the first relatively deep local minimum two particles ap- 
proaching from far away would encounter. It is therefore 
a natural self-assembled configuration for the dimer. Wc 
then pull each along the helical axis with force /. 

At all force levels studied here, the moving dimer ro- 
tates about its centre of mass while d remains perpen- 
dicular to the helical axis (inset Fig. 2e). The behaviour 
at low force (/ < 0.025) is illustrated in Fig. 2(a-d). 
We quantify the rotation by measuring the angle, (j>, be- 
tween d and the x-axis and the angle, 9, between d and 
n(z) (Fig. lc), as a function of time. After an initial 
transient, a smooth rotation is observed (Fig. 2e; open 
symbols), but with separation d that lags behind n(z ) by 
a constant phase angle 9{t) (Fig. 2f, open symbols) [171 ]. 
We attribute this constant lag to a balance between the 
viscous drag opposing the rotation of the pair in the x — y 
plane, and the force arising from the angular variation in 
the rotating interparticle potential. 

At higher force (/ > 0.0275), the behaviour is man- 
ifestly different. The angle between d and x, shown in 
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FIG. 2: [Color online.] (a-d) Snapshots of a steadily rotating 
dimer over a translation of half a pitch (a rotation of 180°). 
The director field is colour coded according to the local direc- 
tor (red along x, blue along y). The time evolution of (e) <j> 
and (f) 6 (defined in Fig. 1C), for the rotor (open symbols) 
and phase slippage (closed symbols) motion. The inset in (e) 
also shows the angle between d and the helical-axis, which 
remains close to 90°. 



Fig. 2e (closed symbols), now increases with time in a se- 
ries of steps. These steps correspond to intervals with and 
without significant rotation of the dimer as it moves along 
the helix (see for a movie and further discussion of 
the dynamics). Correspondingly, the phase lag 0(t) fails 
to reach a steady value (Fig. 2f ; closed symbols) , but in- 
creases indefinitely. We refer to this as "phase slippage" . 
There is a clear transition between the smooth "rotor" 
regime and the phase slippage regime. 

Fig. 3a quantifies the dependence of the average rota- 
tional velocity of the dimer, Q, on the applied force for 
simulations at a range of force values. It can be seen 
that fl is proportional to / in the low force rotor regime, 
while it decreases in the high force phase slippage regime. 
Fig. 3b shows the dependence of the average speed of the 
dimer along the helical axis with force, again showing a 
transition. However, it is difficult to relate these transi- 
tions to the equilibrium potential which is strictly valid 
only at Er = 0. Although the largest Eriksen number 
remains low (Er ss 0.026), the potential is likely to be 
affected by dynamical effects including the local bending 
of the cholcstcric layers (visible in Fig. 2a-d [HI). The 
local bending of the layers for Er <§; 1, was also observed 



for a single colloid moving along the helical axis Q. 

To understand the dynamical transition better, we 
write down a set of phenomcnological equations for the 
evolution of the dimer position, z(t), and for </>(£), which 
gives the direction of d in the lab frame (Fig. lc). We 
assume that the basic features of the cholesteric ordering 
may be captured by an effective angular potential, also 
periodic in z: V((j),z). Our equations read as follows, 
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~dt 
dz 
~dt 



1 dV(<t>,z) 
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-A cos \B(6- 
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Here, A (units of energy) and B (dimensionless) are pos- 
itive constants, / is the external forcing, while 7^ and 
7 Z are relaxational constants related to the rotational 
and translational friction of the dimer, and whose exact 
values we will not need. Eqs. [3] may be solved by an 
ansatz suggested by the behaviour in Fig. 2f. We write 
4>(t) = fit + 0(f), where the director-dimer angle follows 
the most favourable orientation apart from a phase lag, 
0(f). This ansatz is a solution provided that / is smaller 
than a critical threshold f c , i.e. in the rotor phase. By 
estimating the average terminal velocity of the dimer as 
f lltp, we find that in the rotor phase the angular ve- 
locity is f2 ~ 27r// (7^) (this is true within statistical 
error using data in Fig. 3), and that the critical force 
is f c = A'y z p/(2w^ < p). Above this threshold, increases 
with time as in Fig. 2f (solid symbols). The asymptotic 
velocity fi lag (= d0(t)/dt) ~ ^P~f! [u| for / -> f c 
provides a useful "order parameter" to characterise the 
rotor-slippage transition. In contrast, the angular veloc- 
ity 51 decreases as 1/ f for large /. 

Importantly, Eqs. [3] capture both the near-critical be- 
haviour of f2iag, and the large / behaviour of Q shown by 
our full LB simulations (see Fig. 3a for an il ~ 1// fit 
and Fig. 3c for a Vli ag ~ \/ f 2 — f% nt )- The near-critical 
behaviour of Oi ag in our transition is similar to that of 
the velocity of a chain of driven particles in a periodic 
potential described by the Frenkel-Kontorova model [|| , 
suggesting that our dimer provides a liquid crystal rep- 
resentative of a wider class of models . Another ana- 
logue is with the synchronisation of two driven oscillators 
described by the Kuramoto model 0, [TTJ] , where synchro- 
nised and unsynchronised states correspond to the rotor 
and slippage regimes respectively. At the same time, we 
note that the physics of our cholesteric dimers is richer 
than that in Eqs. [3] Whereas the symmetry of the prob- 
lem suggests that there should always be a very low force 
regime in which the dimer behaves as a rotor, the exis- 
tence of the phase slippage regime at higher forces de- 
pends on the form of the effective pair potential. The 
required conditions hold for tangential anchoring but, 
according to our preliminary studies is violated in the 
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FIG. 3: Average rotational velocity (fl) (a) and sedimentation 
velocity (b) of the colloidal dimer as a function of the forcing. 
The non-linear fit in (a) is Q = c/f, with c > a constant, (c) 
shows $lia g as a function of / together with a fit to d\J f 2 — /f , 
with d > a constant and f c « 0.0267 ± 0.0005. 



normal anchoring case, for which we have not observed a 
phase slip regime. 

Our study may be viewed as a generalisation of a clas- 
sical problem: the sedimentation of two spheres in a vis- 
cous fluid. Intriguingly, Fig. 3b shows that the velocity- 
force curves are not linear, even in the Stokes limit of ef- 
fectively zero Reynolds number. Rather, they appear to 
have different slopes (i.e. effective viscosities) in the ro- 
tor and phase slippage regimes. This is different to what 
occurs for a single sedimcnting particle in a cholcsteric, 
which leads to a linear velocity-force relation in the force 
range simulated here This biphasic force- velocity 

curve is due to the dynamic transition we discussed, and 
has no counterpart in classical sedimentation, in cither 
Newtonian or Maxwell fluids. In the Newtonian case, 
sinking side by side speeds up the particles, by up to 



a factor of 2 |18[, which is not true in our rotor phase 



(where we find that a dimer sediments slower than a sin- 
gle particle). In a Maxwell fluid, the repulsive or attrac- 
tive interaction between two spheres sedimenting side-by- 
side is controlled by normal stresses [l9| . In our case, we 
find a novel velocity-dependent torque and a dramatic 
dependence on the nature of the anchoring. 

In conclusion, we have studied the equilibrium and dy- 
namic properties of two colloidal spheres in a cholesteric 
liquid crystal. We have seen that chirality leads to a ma- 
jor change in the effective potential felt by the pair, with 
respect to the nematic limit. The elastic forces we find 
lead to the stabilisation of a dimer at an angle of either 0° 
or 90° degrees, as opposed to the 30° found in nematics. 
These results, alongside those of Ref. Q, suggest that it 
would be instructive to repeat the experiments performed 
in Ref. Q with a cholesteric liquid crystal. Based on our 
results, one may also speculate that variations in particle 
size (or cholesteric pitch) can affect the local free energy 
landscape of colloidal suspensions in liquid crystals, and 
potentially drive the self-assembly of different structures. 

Our main result is that, when subjected to an external 
force, the dimer rotates, either smoothly as a corkscrew 
or intermittently, with phase slippage. This transition oc- 



curs as the forcing exceeds a critical threshold: its value 
may be estimated via a simple theory considering the 
interplay between a spatially periodic angular potential 
and an external driving. The existence of the phase slip- 
page regime is however highly non-trivial and relies on a 
delicate balance in the equilibrium and dynamic proper- 
ties of our dimers: for example we have not observed it 
for dimers with normal anchoring. An interesting possi- 
bility for future research would be to study the effects of 
external electric field applied to our colloidal dimer, as 
done with platelets [23] and nematic colloids (2l[, where 
unusual dynamics was triggered by the field. 
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Supplementary information for: Self-assembly and nonlinear dynamics 
of dimeric colloidal rotors in cholesterics 



COMPUTATIONAL MODEL AND DETAILS 



To describe the thermodynamics of the chiral liquid crystal host, we employ a Landau-de Gennes free energy J 7 , 
whose density / may be expressed in terms of a traceless and symmetric tensor order parameter Q [l[ as 

/ = t(! - i)Ql P - ^QopQmQic + ^fiQl?) 2 

+ ^fopQafi) + "J (ta-fSV-fQsp + 2q Q a p) . (SI) 

Here, Aq sets the energy scale, K is an elastic constant, go = 2ir/p where p is the cholcsteric pitch, and 7 is a 
tcmpcraturc-like control parameter governing proximity to the isotropic-cholesteric transition. Greek indices denote 
Cartesian components and summation over repeated indices is implied; e aj s is the permutation tensor. The physics 
of our cholesteric dimcrs is determined by several dimensionlcss parameters. The chirality, 



V A) 7 



and the reduced temperature, 



t= 27(1-7/3) 

7 

determine the equilibrium phase of the chiral nematic fluid Q. We used the following parameters (in simulation 
units): A = 1.0, K = W ~ 0.065, £ = 0.7, 7 = 3.0, p = 32, q = 1/2 and T = 0.5. Note that in our units the density 
is equal to unity, and so is the rotational viscosity, 71 = Our parameters yield r = 0, re = 0.3. We used a periodic 
cubic simulation box of volume V — 192 3 and particles of radius R = 7.25. Note that the exact forms of Eqs. £[2] and 
331 follows the original work on Landau theory of blue phases by Grebel et al. The chirality, re, differs by a factor 
of 2, from that used by Wright and Mermin Q as explained in the appendix A of For more detailed treatment 
see eg. 

We considered a degenarate planar anchoring of the director n at the particle surface. Following Fournier and 
Galatola , we constructed an orientational order tensor, Q°, planar to the surface normal v as follows: 

Qap = (^7 - v a v 1 )Q 1 s{5si3 - vsvp) - \q8 a p, (S4) 

where 

Q a p = Qal3 + \q5 a f3- (S5) 



CALCULATION OF THE ELASTIC PAIR POTENTIAL 



The calculation of the clastic potential (Fig. la in the main text) was carried out by minimising the clastic free 
energy of equation [T] by solving the Beris-Edwards equations [l| (Eq. (1) in the main text) in the absence of fluid 
flow. This was done for pairs of particles held static at separation d. A first particle is placed with its centre at the 
origin with the helical axis in the z-direction. The far-field nematic director (n(z)) at z = lies in the x— direction. A 
second particle is placed at centre-centre separation vector d from the first. By repeating the simulations for different 
d an effective elastic two body potential can be mapped out. 

We use a discrete representation of the colloidal particle surface in the lattice, which is a standard method in lattice 
Boltzmann models 0, S] ■ This has a techical issue that the volume occupied of a particle can vary depending on 
where exactly the centre of the particle r is located with respect to the lattice. This can have a significant effect when 
considering the elastic interaction energy of colloidal pairs embedded in liquid crystal host. To evaluate the error 
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associated with the discretisation, we carried out a series of simulations where we generated the particle coordinates 
for the colloid pair as, 

ri = r + u, (S6) 

and 

r 2 = n + d. (S7) 

Here ro denotes the origin and the components of u are uniformly distributed random numbers u a £ [0,1], giving the 
position of the centre of the particle relative to the lattice. 

We carried out 8 independent scans consisting of 18 different separations along each of the three coordinate axis 
x, y, z. Here it can be noted that for a constant u, the total volume occupied by the particles (V) is constant for all 
the separations considered along the coordinate axes. 

For the angular scan, we used 19 different angles, 9 £ [0, 180°]. Here, the situation is manifestly different: V varies 
with angle and is symmetric around 9 = 90° only for u a = or 1. This gives a further uncertainty on the angular 
scan in the x — y plane at (d — 2R)/p « 0.05 shown in inset of Fig. la in the main text. From physical arguments 
the clastic interaction energy in the chosen configuration, should be symmetric around 9 = 90°. We confirmed this 
by simulation with u a = 0. However, for u a ^ 0, neither our individual angular scans nor the averaged curve were 
exactly symmetric. So, we further averaged the samples around 9 = 90°, using 16 independent samples for 9 £ [0, 90° [. 

Fig. la in the main text shows the averaged interparticle elastic potential with descriptive error bars of 161fcsT, 
161/sbT and 156/cbT for x, y and z, respectively, whereas for the averaged angular scan all error bars are shown. For 
both cases, these are estimated as the standard deviation of the mean. 



DYNAMICS OF THE DIMER 



To gain some deeper understanding of the dimer dynamics, we replot the time evolution of <fi(t) and 9{t) (Fig. 2e 
and 2f in the main text) against the scaled distance along the helical axis z/p, shown in Fig. £fl] From the 4>(z/p) 
curve the rotor motion (Fig. Sla open symbols) is clear: after an initial transient, the colloid rotates following the 
constant twist of the cholesteric host. 

In the phase slippage regime, the colloidal dimer sediments along the helical axis with non-steady velocity: the 
motion along z alternates between spells of slower and faster sedimentation. The motion is faster when the dimer 
sinks without rotating in the x — y plane. This can be observed by comparing the time evolution of the angle 4>{t) and 
phase lag 9(t) (solid symbols in Fig. 2e and 2f in the main text, respectively) with <fi(z/p) and 9(zjp) (Fig. SI, solid 
symbols). In Fig. SI, intervals of constant <fi denote sedimentation without rotation, whereas intervals of constant 9 
denote rotation locked to the cholesteric helix. By comparing <f)(t) (Fig. 2c) with <f>(z/p) (Fig. Sla), one may observe 
that intervals of rising 4>, corresponding to a roto-translation of the dimer, last longer in time than the constant (f>, 
but covers a similar distance in space - therefore when rotating as a corkscrew the dimer takes longer to move the 
same distance, hence is slower. 

During the sedimentation in the phase slippage regime, the dimer goes through a repulsive barrier in the elastic 
interaction potential at phase lag intervals of 9{t) = 180°. This manifests itself in a repulsion between particles leading 
to a sharp increase of the surface-to-surface separation (r,j — 2R), as shown in Fig. S2 (solid line). For the rotor 
motion, after the initial transient period, the distance between the colloid pair stays constant (dashed line in Fig. S2). 



THE ROTOR-PHASE SLIPPAGE TRANSITION 



Finally, we comment in more detail on theories related to our dynamical transition between the rotor and the phase 
slip regimes. Firstly, our simplified Eqs. (3) may be seen as a relative of the Frenkel-Kontorova model, and of related 
models for the depinning of charge density waves in superconductors. The original Frenkel-Kontorova model considers 
the dynamics of a chain of forced coupled particles in a periodic potential (more realistic models for charge density 
wave depinning consider impurities hence an additional random potential). The key result of the simplest model 
is that there is a transition between _p_inned and moving particles. Just above the critical force, f c , the asymptotic 
velocity of the particles, v, scales as 



v ~ v//2 Pc oc vT^Tc (S8) 
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FIG. SI: Replot of the data in Fig. 2e and 2f in the main text, against z/p. For the rotor motion (open symbols) only every 
5th point is plotted for clarity. 
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FIG. S2: Time evolution of the surface-to-surface gap, nj — 2R, for the rotor (dashed line) and slip motion (solid line), 
respectively. 



where / is the forcing, and 7 is an effective friction. The velocity therefore may be seen as an order parameter for the 
dynamical depinning transition, which is a continuous transition in this framework. While in Eqs. (3) in the main 
text there is a coupling between rotatory and translational motion, the analogue order parameter may be taken as 
the asymptotic time derivative of the lag angle between the dimer orientation and the local director field, 6(i), which 
we call f2ia g - Via a numerical solution of Eqs. (3) we find that f2i ag ~ \/ f 2 — fc V f ~ fc, for / > f c , further 
validating the qualitative similarity hinted at above. In this analogy, our rotor phase corresponds to the pinned state, 
and the phase slip phase corresponds to the depinned one. The critical behaviour appears to be the same, provided 
the appropriate order parameters are chosen. 

Secondly, our model bears some similarities with the Kuramoto model [loj . which is a paradigm to study the 
synchronisation of a system of oscillators. In the particularly simple case of two coupled oscillators, the equations 
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FIG. S3: Plot of fii a g as a function of /, for a numerical solution of Eqs. (3) in the main text (7^, 
simulation units). The solid line is a fit to yjp - /|, leading to f c ~ 6.44. 



1, go = 27r in 



describing the dynamics of the Kuramoto model are 

dffi 



dt 
dt 



= Wl + K/2 [sin (0 2 - 
= u>2 + K/2{sin(6 1 -9 2 )} 



(S9) 
(S10) 



where £q and 02 describe the oscillators, uj\ and w 2 are the forcing for the two oscillators, and K > is the coupling 
constant favouring synchronisation {61 = 62). The relevant dynamical variable is ip = d\ — 6 2l which obeys the 
following evolution equation: 



~dt 



= Au - K [sin (V>)] 



(Sll) 



where Aw = lo\ — uj2- If we define an (asymptotic) angular velocity of the phase difference between the oscillators as 
E = lim^oo then S may be estimated as 2tt divided by the time needed for ip to go from ip to ipo + 27r (this is 
independent of V'o)! namely 



2tt 
At 



2tt 



r 

Jo 



2tv 



dip 

Au-Ksin(tf) 



= v/(Acj 2 - K 2 ) 



(S12) 
(S13) 

(S14) 



where the final integral can be done with the method of residues, by changing variables to the complex number 



e 1 ^ and noting that sin(-0) 



2i 



As in our simplified model, we find that the relevant "order parameter" 



S, behaves as J Aw 2 — AcJJ, with Au c = K, close to the transition. 

Finally, the full lattice Boltzmann simulations of our cholesteric dimers lead to a similar scaling for f2i ag (see Fig. 
3c in the main text). The critical force estimated through the fit is f c sa 0.0267 ± 0.0005, which is in agreement with 
f c G]0.025, 0.0275[ observed in simulations. This shows that our minimal theory captures the physics of the dynamical 
phase transition between the rotor and slip phases, and also suggests that this is in the same universality class of 
oscillator synchronisation and of the depinning of charge density waves in a perfectly periodic medium. 

Note that, although the near-critical behaviour of our "order parameter" is given by ~ \J j — f c , both our simulation 
data and the numerics corresponding to Eqs. (3) in the main text are better fitted by a function proportional to 
\J f 2 — fc as this better captures the behaviour of J7i ag further away from / c . 
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